Unitarity analysis of a non-Abelian gauge invariant 

action with a mass 



In previous work done by us and coworkers, we have been able to construct a local, 
non-Abelian gauge invariant action with a mass parameter, based on the nonlocal gauge 
invariant mass dimension two operator F^,j{D'^)^^F^,^. The renormalizability of the re- 
sulting action was proven to all orders of perturbation theory, in the class of linear co- 
variant gauges. We also discussed the perturbative equivalence of the model with ordi- 
nary massless Yang-Mills gauge theories when the mass is identically zero. Furthermore, 
we pointed out the existence of a BRST symmetry with corresponding nilpotent charge. 
In this paper, we study the issue of unitarity of this massive gauge model. Firstly, we 
provide a short review how to discuss the unitarity making use of the BRST charge. Af- 
terwards we make a detailed study of the most general version of our action, and we 
come to the conclusion that the model is not unitary, as we are unable to remove all the 
negative norm states from the physical spectrum in a consistent way. 



D. Dudai; N. VandersickelJ H. Verschelde* 



Ghent University 
Department of Mathematical Physics and Astronomy 
Krijgslaan 281-S9, B-9000 Gent, Belgium 



Abstract 



1 Introduction 



In the two previous papers IHIIII, the following action was constructed 
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The bosonic fields i?^^, its conjugate and the fermionic (ghost) fields and are 
antisymmetric in their Lorentz indices and belong to the adjoint representation. A"^'^'^ is a 
gauge invariant quartic tensor coupling, subject to a generalized Jacobi identity ||3| 

jman ymbcd _|_ ^mbnyamcd _j_ jmcn yabmd _|_ jmdn yabcra g (14) 

and to the symmetry constraints 

yabcd ycdab 

^abcd ^ ybacd ^ (j g) 

while Ai and A3 are mass coupling^. 

To avoid confusion, let us mention here that we shall work in Minkowski space throughout 
this paper, since we plan to come to the canonical quantization. In (H |2l, the action was 
treated in Euclidean space. 

The classical part of the action, Sd, enjoys a non-Abelian gauge invariance generated by 

















= gr'^'JGl, 




= gr'^u^'G^ 



(1.6) 

with tj" parametrizing an arbitrary infinitesimal SU{N) gauge transformation. 

Quite obviously, the gauge model (|1.1[) did not come out of thin air. Our original motivation 
was based on the quest for a dynamical mass generation mechanism in gauge theories. We 
do not plan to give a complete overview of this issue, but let us mention that this has been a 
research topic since long, see e.g. Q for a seminal work on this. 

More recently, work appeared in which a dynamical gluon mass was introduced phenomeno- 
logically based on the QCD sum rules ||5J. Such a mass can account for ^ power corrections 
in certain physical correlators ||5l[6l[Z|- A natural question arising is where this mass scale 
would originate from? The authors of (6l[7l invoked the condensation of the operator 

^Ln = (VT)-' ^ min^^ / d'x (^^)^ (1.7) 

since it is gauge invariant due to the minimization along the gauge orbilH. As it is well 
known, a local gauge invariant dimension two operator does not exist in Yang-Mills gauge 
theories. The nonlocality of (|1.7|) is best seen when it is expressed as a series in Euclidean 
space Cni 



— / fJ^T 

^min ~ 2VT 



+ ..(1.8) 



^In comparison with 1 1 2 1, we changed the sign of m, A3, Ai and A"'"^'' to avoid a number of minus signs. 
^One should however be aware of the problem of gauge (Gribov) ambiguities 1 8 9 1 for determining the global 



mmimum. 
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which contains the inverse Laplacian ^ several times. This is a nonlocal operator, as it can 
be immediately inferred from its formal expression in d dimensions through 






fiy) 



(1.9) 



27r2(d- 2) 



x-y 



d-2 ■ 



X 



All efforts so far were concentrated on the Landau gauge d^A^j^ = 0. The preference for 
this particular gauge is obvious since the nonlocal expression (|1.8|) reduces to an (integrated) 
local operator, more precisely 



In the case of a local operator like A^, the Operator Product Expansion (OPE), viz. short dis- 
tance expansion, becomes applicable, and consequently a measurement of the soft (infrared) 
part (^^)opE becomes possible. Such an approach was followed in e.g. fTT\ by analyzing the 
appearance of ^ power corrections in (gauge variant) quantities like the gluon propagator 
or the strong coupling constant, defined in a particular way, from lattice simulations. Let 
us mention that already two decades ago attention was paid to {A'j^)opE when the OPE was 
applied to the propagators [l2J. This condensate (A^)ope can also be related to an effective 
gluon mass, see e.g. 1)13 J . 

A more direct approach to a determination of {A"^) in the Landau gauge was presented in 
flUllS^ . In |[14|. a meaningful effective potential for the condensation of the local composite 
operator ^4^ was constructed, giving evidence of {A"^) / 0, and as a consequence a nonvan- 
ishing gluon mass of a few hundred MeV was found. The renormalizability of this technique 
was proven to all orders of perturbation theory in |ll6| . 

Effective gluon masses have found application in phenomonological studies like |fT7l[T8l[T9| . 
Also lattice simulations of the gluon propagator revealed the need for massive parameters, 
when the obtained form factors are fitted by means of functional forms [20, 21112211231. Other 
approaches to dynamical gluon masses are e.g. Il24l |25| . The estimates of these mass pa- 
rameters are grosso modo all in the same ballpark, ranging from a few hundred MeV up to 



It is perhaps important to spend a few words at clearing up a common misconception. The 
concept of a dynamically generated effective gluon mass does not necessarily entail that we 
are considering massive gauge bosons that are belonging to the physical spectrum, i.e. that 
are observable particles. At low energies, perturbative QCD expressed in terms of gluons 
and quarks completely fails, and the effective degrees of freedom become the hadrons. The 
phenomena we are interested in, in casu the study of the condensates and ensuing dynam- 
ical mass generation, occur in a energy window located in between perturbative QCD and 
the confined region. Perturbation theory still has its validity there, but it gets corrected by 
nonperturbative effects like condensates. Due to the lack of an explicit knowledge of the 
correct physical degrees of freedom (the hadrons), we continue to use the gluons as effective 
degrees of freedom, although we are already out of the energy regime where these might be 
considered as asymptotic observables. If we cross from high to low energies, the originally 
massless and physical gluons will not stante pede become confined at the confinement scale, 
but rather they will behave as a kind of massive quasi particles before getting confined, and 
this happens at scales that are phenomenologically relevant. This also means that unitarity 




(1.10) 



1.2 GeV. 
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Gauge 


Operator 


linear covariant 

Curci-Ferrari 
maximal Abelian 
nonlinear class 


i Aa Aa 

\A^X + orc^^^ 

1 AO. Aa 



Table 1: Gauges and their renormalizable dimension two operator 



in terms of the gluons is not required or even desired. One expects that quasi particles do 
have a finite lifetime and cannot be observed as asymptotically free particles. 

We have already explained the preferred role of the Landau gauge, since in that case a gauge 
invariant meaning can be assigned to (A^). Obviously, since we are working in a gauge 
theory, the condensates influencing physical quantities should be at least gauge invariant. 
Therefore, it would be nice to have a dimension 2 condensate that could also be treated in 
other gauges. As the operator A^j^^ remains nonlocal, it falls beyond the applicability of the 
OPE. It is also unclear how e.g. renormalizability or an effective potential approach could 
be established for nonlocal operators. In most covariant gauges, we and collaborators have 
discussed that other dimension two, renormalizable and local operators exist. We showed 
that these operators condense and give rise to a dynamical gluon mass, see Table 1 and 
|l26l|2Zl|28l|^|30H3TJ|32l|33l. Quite recently, it has also been shown that a class of nonlin- 
ear covariant gauges enjoys the fact that is multiplicatively renormalizable ||34| . In the 

maximal Abelian gauge, it was found that only the off-diagonal gluons A^^ acquire a dy- 
namical mass [31 J, a fact qualitatively consistent with the lattice results from p2ll23] . Let us 
also mention that we have been able to make some cormection between the various gauges 
and their dimension two operators by constructing renormalizable interpolating gauges and 
operators I3TII35| . These can be used to obtain a formal result on the gauge parameter in- 
dependence of the nonperturbative vacuum energy due to the condensation, which is lower 
than the perturbative (zero) vacuum energy Il30l . 

A certain disadvantage of the research so far is the explicit gauge dependence of the used 
operator. We started looking for a gauge invariant dimension two operator, which a fortiori 
needs to be nonlocal. We would like to develop a consistent (calculational) framework, hence 
we are almost forced to look for an operator that can be localized by introducing a suitable 
set of extra fields. From this perspective, A^:^^ seems to be rather inadequate as it is a infinite 
series of nonlocal terms. A perhaps more appealing operator is HI 

1 



VT 



ab 

(1-11) 



This operator found already use in the study of a dynamical mass generation in 3-dimensional 
gauge theories [36J. When we add the operator O to the Yang-Mills action via 

t2 



Sym- — J d'xF^,[{DY ] F'^,, (1.12) 



we can localize it to 



Sym + I d'x {B - B)l^ f;, + i [B^DfD'fB^, - G%DfD'^^G%] 



,(1.13) 
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at the cost of introducing a set of extra fields |fT|. 

The action (|1.13|) as it stands is however not renormalizable, but we and collaborators have 
shown that the generalized version is renormalizable to all orders in the class of lin- 
ear covariant gauges, implemented through Sgf, in We have also calculated several 
renormalization group functions to two loop order, confirming the renormalizability at the 
practical level. Various consistency checks were at our disposal in order to establish the reli- 
ability of these results, e.g. the gauge parameter independence of the anomalous dimension 
of gauge invariant quantities like g^, \°-^'^'^ or m, or the equality of others, in accordance with 
the output of the Ward identities in [IJ. We refer the reader to [Tl|2] for all details concern- 
ing the localization procedure or renormalizability analysis, as well as the need for the extra 
couplings. 

Furthermore, we have proven in |j2l the perturbative equivalence of the model with or- 
dinary Yang-Mills theory in the case that m = 0. We notice that this is a nontrivial statement 
due to the presence of the quartic interaction ~ X"-^(^d- in the extra fields. It has an interesting 
corollary: because we employ a massless renormalization scheme, in casu MS, we can set the 
mass m equal to zero to determine the renormalization group functions of e.g. the coupling 
constant , the gauge parameter ^ or original Yang-Mills fields. Since both theories are per- 
turbatively equivalent for m = 0, the already mentioned renormalization group functions 
must be identical. This has indeed been confirmed by the explicit results of [1 , 2J. In particu- 
lar, our model is thus asymptotically free at high energies, with or without a mass. At lower 
energies, nonperturbative effects can set in, completely analogous to the Yang-Mills case. 

Summarizing, we have thus found a classically gauge invariant action, which at the quan- 
tum level can be renormalized to all orders in at least the class of linear covariant gauges, 
and as a bonus it is perturbatively equivalent with ordinary Yang-Mills gauge theories for 
vanishing mass. We can now ask ourselves two questions: 

1. If we treat the mass m as a given classical input, can we consider our model as a candi- 
date for a gauge theory with massive excitations? Therefore, we should prove that the 
theory is unitary, containing massive particles in a suitably defined asymptotic phys- 
ical subspace. The particles correspond to the elementary excitations of the original 
fields. As it is well known, proving the unitarity of gauge theories is not a trivial job. 
A well known proof in the case of Yang-Mills theories based on the BRST symmetry 
EZIIMI, is given in l.39il. 

2. If we do not want to treat our model as one with a given classical mass m, can we 
dynamically generate it in a selfconsistent way in this case? Said otherwise, can we de- 
velop a method to find a reasonable gap equation for this mass? At high energies, the 
model is massless and the same as Yang-Mills theory, but it might develop a dynamical 
mass scale at lower energies, without spoiling the gauge invariance. We cannot add 
mass terms to the Yang-Mills action without spoiling the gauge invariance or renor- 
malizability, but we can add mass terms to our model. Just as for Yang-Mills theories, 
we expect our model to be confining at lower energies. As we have already mentioned, 
the original fields can develop a behaviour different from the one expected from per- 
turbation theory in the energy regime in between confinement and the perturbatively 
accessible high energy region. For example, a gauge invariant mass parameter could be 
dynamically generated, thereby modifying the propagators in a nonperturbative fash- 
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ion, and this without the need that these describe asymptotically observable physical 
particles. 

In this paper, we shall provide an answer to the first question by studying the massive gauge 
model (|1.H) . More precisely, we shall quantize the model canonically to have a clear particle 
interpretation of the quantum fields, and we shall find out whether it is possible to define 
a physical subspace of states endowed with a positive norm. Naively, one might expect 
the model to be unitary, because the action (|1.H) enjoys a BRST symmetry, generated by the 
nilpotent operator s, 











si?" 


















= }f , 


sh" 


= 0, 




= 0. 



It should not come as a surprise that we shall heavily rely on this BRST symmetry to dis- 
cuss the unitarity of the model. The paper is organized as follows. In section 2, we review 
how a sensible physical subspace can be defined by using the free BRST charge |l40ll4]l . As 
a warming up exercise, we apply the results of section 2 to the well known canonical quan- 
tization of Yang-Mills gauge theories in section 3, before turning to the explicit quantization 
of the gauge model il.l} in section 4. In section 5, we discuss the presence of some extra 
symmetries which allow to reduce the physical subspace further. Section 6 is devoted to the 
(free) classical equations of motion and the Fourier decomposition of the solutions. We shall 
encounter the problem of "multipoles", since the free equations of motion couple different 
fields to each other. This gives rise to higher derivative decoupled equations of motion. We 
also pay attention to the BRST charge and Hamiltonian. In section 7, we discuss how to 
derive the commutation relations between the creation and annihilation operators, without 
using the brackets between the fields and their conjugate momenta, which we want to avoid, 
since not all the Fourier components of the fields and momenta are independent. Once this 
is done, we come to the conclusion that the massive gauge model (ll.lj) is not unitary, as we 
end up with negative norm modes in the physical subspace. We are unaware of any step to 
further reduce this subspace in a consistent zraijl to remove these unwanted modes. We end 
with some conclusions in section 8. 

■^That means compatible with the interactions of the model. 
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2 A constructive approach to the question of unitarity in gauge the- 
ories 



In this section, we shall review how we can construct the action S for an interacting gauge 
theory, if we have a free theory at our disposal, together with a nilpotent symmetry gen- 
erator So, so that sqSq = 0. The content of this section is mainly based on [40, 41], although 
here and there we adapted the proofs. We shall only be concerned with the non-reducible 
case in this paper. 

Let us thus start from the free action Sq. This action contains a set of fields, appearing 
quadratically. It is given that enjoys a BRST symmetry sq, with corresponding nilpo- 
tent charge Qq. As a standard example, we can consider the free part of a gauge theory 
in a particular gauge, with its corresponding gauge fixing part. For example, in the linear 
covariant gauge we have 



Sn = d X 



(2.1) 



The free BRST symmetry is generated by 



SQC^ = 0, 

so'f = b\ 
sob"" = 0, 



si = 0. (2.2) 

We mention that also the free "ghost part" has to be included in 5*0. We may define the 
ghost charge Qgh- In (40ll4l1l , the ghost charge is not used. We may use it anyhow in the 
definition of the physical subspace. However, this requirement is a bit redundant. A BRST 
cohomological analysis (see later) will eventually learn that a physical state counts neither 
ghosts nor anti-ghosts in the case of Yang-Mills gauge theories. 

Unitarity means that we start from a physical state space Wphys, which is a subspace of the 
total Hilbert state space H. H^^ys should of course be endowed with a positive norm in order 
to have a sensible probabilistic interpretation of the quantum theory. If we let the states 
of Wphys interact, we must end up again in the (same) subspace Wphys- Nonphysical states, 
which can have negative norms, may contribute to the 5-matrix in internal processes, but 
they cannot appear in the observable sector (the "out"-space), unless perhaps in zero norm 
combinations. 

Consequently, two questions need to be answered: 

1. How do we define the physical subspace H^^yP. 

2. Do the states in the physical subspace H^hys possess a positive norm? 

Let us first explain how we define our physical subspace, starting from the free action. A 
state IV'p) is called physical if 

IV'p) eWphys^Qo|V'p> = 0, \i;p)^\...) + Qo\...). (2.3) 
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Physical states are thus defined from the "free" BRST charge Qq. Since Qo is supposed to 
be nilpotent, states of the form Qol • • •) are trivially annihilated by Qo- We notice that these 
have zero nor We can identify them with the trivial state, more mathematically speaking 
this amounts to consider the Qo cohomology. In the usual terminology, we define Qo -closed 
and Qo-exact states by 

iVp) is Qo-closed ^ Qo|V^p) = 4^ iV'p) G KerQo , (2.4) 
iV'p) is Qo-exact ^ \^p) = Qo|(/>) ^ \ipp) G ImQo . (2.5) 

Since Qg = 0, every exact state is trivially closed, meaning that ImQo C KerQo. Hence, we 
can reexpress the condition (|2.3|) as 

KerQo 

IV'p) e Wphys ^ iV'p) G cohomQo = ^ ^ . (2.6) 

iiriyo 

For the moment, we leave open the (key) question whether these states iV'p) have a positive 
norm. 

The next problem is whether we can construct an action S compatible with unitarity? Start- 
ing from the free action So, we can complement it order by order with terms in the coupling 
constant(s), so that 

S = So + Si + S2 + ... . (2.7) 

The question becomes how to determine the interaction terms 5*1, 5*2, such that S describes 
a unitary model? More precisely, having defined a physical subspace by means of (|2.6|), we 
would like to construct the action 5 such that the subspace defined by (|2.6|) is maintained 
under time evolution. In the operator language, we must therefore require that the time 
evolution operator S, given by 



r 



(2.8) 



with T the usual time-ordering operation, commutes with the operator Qo. Then clearly the 
5-matrix will be unitary, as states evolved w.r.t. S will asymptotically again belong to the 
(same) physical subspace Wphys- 

In order to solve the previous requirement, we prefer to work in the path integral language 
rather than in the operator language. Let us thus rephrase the previous requirement in the 
path integral language. From the LSZ reduction formulae, see [42J and [43J for the original 
paper, we know that the 5-matrix element^ are determined by the (connected) amputated 
n-point Green functions, put on-shell. In a rough notation, we can write 

(fc'l, . . . , k'JS\ki, . . . , 4) ~ (e|T [(t>,r.(t>n] |e) , (2.9) 

where ~ symbolizes all the necessary prefactors, putting it on-shell, amputating and Fourier 
transforming to momentum space, (pm and ^„ are certain functionals of the fields, leading to 
a (m + n)-point function. 

Starting from a generic physical state \ki,. . . ,kn) with the property 

Qo|^i,...,4) =0, (2.10) 



*The BRST charge can be chosen to be Hermitian. 

*We can restrict ourselves to the connected 5-matrix elements. 
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we are wondering which condition will assure that 

{k[, . . . ,k'^\QQS\ki, ...,kn) = {k'l, . . . ,k'^\SQo\ki, . . . = 0, (2.11) 

where the last equality follows from l|2.10|l . As it is well known, we can express the T-ordered 
product with the path integral, so that we finci 

{k[,...,k'JQoS\ki,...,kn)-{e\T[{Qo(Prn)Mm = J dX {Bo^m) cj^ne'^ , (2.12) 
where X represents all the fields. Now, we can write 

J dX (Bo^rn) ^ne'^ = J dXBo {(l)m<Pn) e'^ , (2.13) 

since Bo4>n = by virtue of (|2.10|) . 
Consider the path integral 

' dX^e'^ , (2.14) 



I' 



for an arbitrary functional <I> of the fields. We perform the transformation of the path integral 
variables 

X = X' + eBoX' , e infinitesimal . (2.15) 

As Bq induces a linear transformation, there is no associated Jacobian, and we find 

dX^e*'^ = j dX''^'e'^' + e j dX'Bo (^$'e^^') . (2.16) 

Dropping the prime again, we find 

j dXBo ($e^^) = . (2.17) 

Taking a look at (|2.17|) , we can be certain that (|2.13|) holds when we impose 

Boe'^ = . (2.18) 

In order to proceed, we notice that it is in principle sufficient that (|2.18|) is fulfilled on-shell 
as the 5-matrix is of course considered on-shell. At the level of the action however, we must 
require that it holds off-shell. Let us introduce a (very) condensed notation for the action of 
Bo 

Bo = A0o^ ■ (2.19) 
Implementing (|2.18|l at lowest order and making it valid off-shell means that 

BoSi = A0o7r7^i = -A0i^. (2.20) 

od) 06 

off-shell 



''We shall use the notation Qo for the charge, eventually written in terms of creation/ annihilation operators, 
while Bo represents the fimctional analog of Qo • 
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We already see here that Bq will get adapted, more precisely we can introduce the modified 
operator B by 

Bo ^ A0o^ -^B = Bo + Acl>,^, (2.21) 

so that 

B{So + Si)=0. (2.22) 
We remind here that all Q's (B's) are Grassmann operators. 
Since Bq is nilpotent, we can act with it on (|2.20|) to find that 

= ^2^1 = ^0 f-A</.i^5o 



Acting with ^ on 
yields 



+ ^'^i^'^ol^^o . (2.23) 



= BoSo = A(l)o—So , (2.24) 



^(A0o)^5o + A0o^5o = O. (2.25) 



Combination of l|2.23|) and (|2.25|l learns 



from which we infer that 

The identity (|2.27|) expresses nothing more than the nilpotency of Q = Qo + Qi/ given by 
(|2.21|) since, at lowest order 

where we used for example the nilpotency of Bq. We dropped the last term as it is of higher 
order. 

We notice that the potential solution of (|2.20|) is apparently restrained by the condition that 
it is invariant under a nilpotent operator Q (B), which reduces to Qo (Bo) in the free limit. 

This construction can be continued at higher order. One proves that the action at order n, 

S = So + Si + ... + Sn, (2.29) 
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is the solution of 

A</.o^5„ + A^i—Sn-l + ... + ^^n-g^So = , (2.30) 
where consistency demands that the BRST operator, 

B = Bo + ... + Bn = A4>o— + ... + A0„— , (2.31) 

is nilpotent at the considered order n, thus 

= (or g2 = 0) . (2.32) 
By construction, the final action (|2.29|) shall be invariant under the BRST symmetry generated 

To make things a bit more comprehensible, let us work out the procedure at second order. 
We hence demand that (|2.7|) is consistent with (|2.18|) , and this extended to the off-shell level, 
meaning that 

iBoS2 - -BoiSiSi) = -iA^2^ . (2.33) 
2 0(p 

The complex unity i in the r.h.s. as well as the ~-notation are merely introduced for later 
convenience. Using (I2.20D , we may rewrite ( I2.33D as 

iAcj)o^ + A(/>i-f 5i = -jA02^ . (2.34) 
ocp ocp ocp 

Next, using Wick's theorem, we can writ^ 

. , dSn ^ , . , dSi — -- dSn 

A^,^5, = zA^,— +.A0,^, (2.35) 

since roughly said, ~ "D"^ x 0, and a "contraction" of this with a (/> from S\ will give rise 
to a iD^^ X V with V a free propagator. All other terms are taken together in A0i. 
Upon taking (|2.34|) and (|2.35|) into account, we come to the conclusion that 

^*«t^^*'t-^*=t ''''' 

in order to have the condition (|2.33D fulfilled. We defined 

A02 = A^2 + A^i • (2.37) 

Analogously at is was proven in (9) to (14), the nilpotency of Qo + Q\ leads to the nilpotency 
of (5o + <3i + Q2 as a consistency requirement. 

Of course, there is no guarantee that the foregoing "bottom top" construction of the complete 
action will end at a finite order. Given that it ends at a finite order, it could still be a very 
cumbersome job to actually get the nilpotent BRST charge and corresponding action. The 
situation becomes much more appealing when we already have at our disposal a complete 



^This operation is understood within the path integral. 
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action, with a nilpotent charge generating a symmetry. If the interaction is switched off by 
setting all coupling constants equal to zero, we obtain the free action, with a free nilpotent 
charge. When the above "bottom top" machinery is unleashed, the complete original action 
and its BRST symmetry generator shall quite evidently be a solution to the iterative proce- 
dure. From this viewpoint, we have a "top bottom" approach to unitarity for actions with 
nilpotent BRST charge, when they are "reduced" to their free counterpart. 



3 Unitarity of Yang-Mills gauge theories using the BRST charge 



We should still provide an answer to question 2, namely do the states that are annihilated by 
the free BRST charge have a positive norm? It is well known that this is the case for Yang- 
Mills gauge theories. For completeness, let us nevertheless repeat the argument. This will 
allow for a comparison with Yang-Mills theories when we start analyzing our generalized 
model. 

We shall base ourselves on 144) for this particular jolfl We opt to work in the Feynman gauge 
for simplicity = 1 in l|2.H) ). Let us first determine the conjugate momenta of all fields. 



vr 



0,a 



77a 

-doc^ , 



so that quantization requires 



A^{x,t),7r'jy,t) 



'A-,{x,t)y'''{y,t) 
{c»(x,t),7r,^(y,i)} = 

other (anti-)commutators trivial . 
We mention that the classical equations of motion are 



iS'^H'^^^x - y) , 
iS'^H^^^x - y) , 



,2 AO, 



Qlia 



0, 



0, 



and that we use the hermiticity assignment 



(3.1) 



(3.2) 



c, 



(3.3) 



(3.4) 



^We shall however use other conventions than those of l44l 
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We propose the following Fourier decomposition^ 

d^k 1 



^■>(-> = /72rf2:r +«S(*)^'" 



m=l 



"("^^ = /(03^(^(^)e"^'"-^^(^K'') • (3-5) 

The polarization vectors ef"^ (fc) form an orthonormal set, ef^^ (^) = <^'"" with e^^^^ (fc) = 

1^ . We shall assume that the particles move along the z-axis, so that = Sij . 

Implementing (|3.2|) , we must require the following (anti-)commutation rules 

[ao{k),ao{q)] = -2{27Tfu;kd^^\k - q) , 
[a^{k),an{q)] = 2{2^fuOk5^^\k - q)5mn , 
{^{k),r,\q)] = -2(27r) - , 

{r/(A:),r/t(g)} = -2{2TTf u:k6^^\k - q) . (3.6) 
For later use, let us already introduce the operator 



{-al{k)ao{k) + alik)as{k) - 7^\k)m - fj^k^k)) , (3.7) 
which counts the unphysical modes. 

We are now ready to express the BRST charge in terms of the creation/ annihilation operators. 
The BRST Noether current is given by 

^ Ff^^'^d^c" -b^df.c", (3.8) 

which leads to the charge 

Qo = J d^x (c"5°6" - h^'d^d') , (3.9) 
where use was made of the classical equation of motion 

After substitution of (|3.5|) in (|3.9|) , the BRST charge is expressed as 

d^k 



Qo 



{27ry 



al{k) - alik)) r]ik) + r^^{k) (ao(A;) - a^ik)')] . (3.11) 



'We suppressed the global color indices. 
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alik) + alik)) r]ik) + r]^{k) (ao(fe) + a3(A:))l , (3.12) 



If we define 

^ _l [ d^k 1 
2 J l2^J^l 

then a little algebra yields 

AA={Qo,7^}. (3.13) 

The fact that the "nonphysical" counting operator J\f is BRST exact is a very powerful result 
[44|. Assume that IV'p) is constrained by 

Qo|^p)=0, (3.14) 
and that it contains n 7^ unphysical modes, i.e. 

AAlV'p) = nlV-p) , (3.15) 

then consequently 

\^p) = -\^p) = - {Qo-R + nQo) = Qo (-nii^p)) , (3.16) 
n n \n J 

meaning that a state |^p) annihilated by the BRST charge and containing nonphysical modes 
is a fortiori BRST exact, and hence it is zero in the physical cohomology. Said otherwise, 
physical states do not contain unphysical modes. The physical subspace Hphys of Yang-Mills 
gauge theories does only contain the 2 transverse polarizations of the gauge field, whereas 
the scalar and longitudinal polarizations cancel with the ghost degrees of freedom. 

In (lot , s different proof was presented of the fact that a state annihilated by a BRST charge 
Qo of the form (|3.1H) containing nonphysical modes, must have zero norm. However, the 
cohomological approach used in e.g. [ (44,1 is somewhat more elegant. 

4 Application to the massive gauge model: preliminary remarks 

Setting the couplings g and X"-^'^'^ equal to zero in (|1.1|), we are considering the quadratic 
action 



So = d X 



+ m^^ {B% - B%f + b'^d.A^^'^ + ^d\- + 
This action enjoys the free BRST symmetry generated by 



p,v — •sqC^j, — SqG^i, — , 



(4.1) 



soAl 


= -d. 


soc" 


= 0, 




= soB'i 


soc^ 


= h\ 


sob" 


= 0, 



(4.2) 
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where clearly 

sl = 0. (4.3) 

We can hence apply the results of section 2 to the action The only thing left to prove 
is that there exist a physical subspace with positive norm. This subspace is certainly anni- 
hilated by the (free) BRST charge, but nothing prevents us from using other available sym- 
metries to further reduce the physical subspace. In the next section, we shall introduce 2 
extra symmetries with nilpotent generator of the complete action We first determine 
the BRST charge in functional form. We shall see that it remains unchanged compared to the 
Yang-Mills case (|3.9|l . The Noether current corresponding to the BRST transformation (|4.2|) 
and action (|4.1|l is given by 

= F^"''''^^c''-b''^^'c''-'-^{B-B)'"''''^^c'', (4.4) 
which leads to the BRST charge 

= /.3.(e.,0,._,.,0,.), ,,,, 

where we invoked the equation of motion 

Q^p,u ^ QU^a ^ 1^^^ _ _ ^^ ^^ 

For what concerns the Faddeev-Popov ghosts c and c, it is immediately seen from the action 
( 14.11 ) that their quantization remains unchanged compared to the Yang-Mills case, given in 
(|3.H) , (|3.5|) and (|3.6|) . Therefore, since (|4.5|l must be time independent as a conserved charge, 
we already infer that Qo will only act nontrivially on massless excitations. This shall be 
confirmed later once we have found the excitations belonging to the 6" -field (see section 6). 



5 A further reduction of the physical subspace 

In the following sections, we shall make use of a cohomological result (451, summarized 
here. 

Doublet theorem Consider a transformation b with the property that 

6v^ = 0, 6v', = 0, ^""-'^ 

with Ui,v^ commuting and u'^,Vi anticommuting quantities. We call {ui,Vi) and {u'-,v'j) 5- 
doublets. 

Then it is a trivial exercise to show that 5 is a nilpotent transformation. Moreover, (uj, Vi) 
and (n-, u •) appear trivially in the (5-cohomology. This can be proven 1.45 J by introducing the 
"counting" operator 
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and the operator 
such that 

V = {6,A} [P,6]=0. (5.4) 

Assuming that 5X = 0, we can expand X in an mutual eigenbasis of the commuting Hermi- 
tian operators V and 6. It is then quite easy to show that 

X = Xo + 6Y, (5.5) 

whereby VXq = 0, i.e. the cohomology of 6 does not depend on (uj, Vi) and (n-, f ■). 

In the paper |l2|, we already noted the equivalence between the action il.l^ for m = and 
conventional Yang-Mills theories, when quantized using the same gauge fixing. The physical 
content of both theories should be the same, in casu the unitarity should be satisfied. This 
can be shown in the following way. The action ^1.1} enjoys a supersymmetry when m = 0, 
generated by the nilpotent transformation [2J 

= 0, 6s%, = Q. 

So, the Doublet theorem applies, and we conclude that the excitations belonging to the extra 
fields will not belong to the physical subspace. 

In the case that m 7^ 0, the supersymmetry \5.6) is broken due to the terms ~ {B — B). In 
a matter of speaking, the symmetry is only broken by terms ~ (i? — B), and not by terms 
~ {B + B). Therefore, we might expect that some trace of the supersymmetry 5s might 
survive after all. We shall first explore this possibility. We decompose the fields G^^, and G^^, 
in their "electric" and "magnetic" part 

u. — 2^ jfc ' . ~ 2^ '^jki (^7\ 
pi,a _ Qoi,a '0''°' ^o*i<^ ^ ' ' 

Consequently, one finds 

- ^l^^'^G^"'^'' = ^f'^'d'^P''^ - hx^'^d^a''^ . (5.8) 
Since G and G always appear in the "product" combination G't^,0'^^G^'^''^ , only the structures 



j30(j and aOa will appear. Having a look at the complete action \\.\} , it is clear that also in 
the interacting theory, these are the only possibly appearing structures. 

It is convenient to also decompose the B and S-fields in their electric and magnetic counter- 
parts 

_ -eJ B-f^, p' - Bj,,, 
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in which case the (quadratic) action becomes 



2 2 2 2 2 



—t,a—t,( 



8 ^ ^ 16 ^ 



(5.10) 



For further convenience, we shall exchange {a, a) and {p, p) for their real and imaginary- 
parts via 



(5.11) 



yielding 



5o 



2 2 2 



— -m' 
4 



As a final step, we introduce the fields 



,±,i,a 



X 



(5.12) 



(5.13) 



to write 

^0 



+ l^+,i,a^2^-,i,a _ ljYi,a^2^i,a ^ l^i,a^2^j,a ^ l^."^2^i,a _ \-^i,aQ2^i,a 

2 2 2 2 2 



We introduce the following transformations 



(5.14) 



(5.15) 
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and 



:,a 



(5.16) 



Clearly, these transformations define a symmetry of the free action (|5.14|l . However, these 
will also generate a symmetry of the full action jl.l} . One notices that only the 3rd, 4th and 
6th term of (|1.2|) will give rise to contributions in x^- Taking a closer look at these terms, it is 
quite easily seen that these contributions will always be of the type 

^+,t,a(jab^-,i,b ^ ■^,a^ab^i,b _ -i,aQab^i,b ^ ^5 -^7^ 

with 



and one can check that {5.V7) vanishes when ,3 is applied to it. 
It is also readily derived that 

{A,,A^} = 0, (5.19) 

while it also holds that 

{A„,s} = = {A^,s} , (5.20) 

since the BRST operator s acting on the new fields can be read off from dl.MI ), (|5.7|) , (|5.11|) 
and (I5l3ll to be 

^a,i ^ |a^.'^,a*''^,/3''^;5^''',x+'''",X"'*'°,^''",A^''''} ■ (5-21) 

Recapitulating, we have found 2 symmetries of the action dl.ll) which are generated by the 
nilpotent generators (|5.15|) and (|5.16|) . Moreover, the fields (a* ", x^'''") ^rid (x a*'") form 
Aq doublets, so we can be assured that these fields decouple from the physical sector by 
applying the Doublet theorem. Moreover, we can equally well deploy the Doublet theorem 
on A^ to also remove [3 and [3 from the physical subspace. The intersection cohom(Aa) n 

cohom(A/3) is thus independent of a*'", a*'", x"*^'*'" and x~'''"- 

Consequently, we can conclude that all degrees of freedom corresponding to the extra ghost 
fields and as well as the extra bosonic degrees of freedom corresponding to + 
B""^^ are decoupled from the physical sector. 

For the fields A^*''* and L*'", corresponding to the degrees of freedom in i?^^, — B"'^^, the 
analysis is less quickly performed, as these are coupled to the original gluon field in 
a nontrivial way. In the next sections, we will have a look at this problem. As we did 
not analyze yet the space annihilated by the free BRST charge, we may expect that certain 
degrees of freedom will also be killed when this subspace is considered. 



ja6 
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6 The classical equations of motion and the Fourier decomposition 
of the fields 



As a next step, we must quantize our model. Before turning to the conjugate momenta and 
quantization rules, it is advisable to have a look at the classical equations of motion in the 
new variables and the corresponding Fourier decompositions of their solutions. Due to the 
mixing between the fields, their Fourier coefficients will not all be independent. 

6.1 Classical equations of motion 

The free classical equation of motions arj^ 



b + d^'A^ 


= 0, 


(6.1) 


d'^Ao - md^N' 


= 0, 


(6.2) 


+ m'^)N' - m(a°A* - d'A^) 


= 0, 


(6.3) 


-{d"^ + m^)n -me'^'^djA^ 


= 0, 


(6.4) 


d'^A' + mdoN' - me'^^djL^ 


= 0, 


(6.5) 



where we set 

= — — m^Ai — — m^Aa . (6-6) 
Since these equations are coupled, the quantization of the fields is not straightforward. 



6.1.1 Intermezzo: multipole fields 

After some manipulation with the equations (|6.H) - (|6.5|) , we derive that 

(a^ + m2 + (S^Ao) = 0, 

52 + m^ + fh?) (diN^) = , 

{d^+m^)idiV) = 0, 

d'^ {d^ + m'^ + m'^)d^{diA') = 0. (6.7) 

These (induced) equations of motion are of higher order in the derivatives. In the literature, 
such fields are known as "multipole" fields, and their quantization is indeed more involved 
than in the well known Klein-Gordon case. The problem of dipole fields also appears in the 
supergravity context Il46| . 

Consider e.g. the following toy model of a higher derivative action ||46| 

S = -h^d'<P, (6.8) 

with classical equation of motion 

d'^d'^cf) = . (6.9) 



'^We shall skip again the global color index. 
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It is not clear how to quantize the model (|6.8|) since higher order derivatives occur. We can 
rewrite (I6.8I) as 



S 



1 



(6.10) 



by introducing an auxiliary field iJj- It is important to notice that the derivates now occur at 
most quadratically, making the action \6.\Q) suitable for canonical quantization. In fact, this 
is the case we are investigating. Our action \5.12) is indeed already at most quadratical in 
the derivatives. 



The equations of motion associated to (|6.10|) are 



32„ 



9> = 0. 

Let us turn to Fourier space. We propose the solution 



(27r)3 
(2vr) 



V'(g,t)e*«^ + h.c., 
3^(g,t)e^«^ + h.c.. 



(6.11) 

(6.12) 
(6.13) 



Plugging \6.12) in (|6.1H) , the following differential equations in the time t = xq arise 

{dl + if)^{q,t) 



m'^'tp{q, t) 
0. 



We can set 



^P{q,t)=a{q)e-"^^\ qo = \q\ , 
hence the equation l|6.14|l is satisfied when 

,2 



m 



t) = 6(gle-*'?o* +-^(1 + 2i\q]tMQ)e- 



-iqot 



Aq 



qo = lol 



(6.14) 
(6.15) 

(6.16) 



solution hom.eq. 



According to |[46ll , terms linear in time t are typical for dipoles. As we shall soon see, terms 
oc t shall also appear in our massive gauge model. One can then continue by quantizing the 
model. The necessary commutation relations are of the type 



a{q),b\k)] = \b{q),a\k)] = -Lsi^) - ^ . 



(6.17) 



Similar, although a little more complicated techniques and relations, will occur when we try 
to solve the classical equations of motions of our model in Fourier space. 
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6.1.2 Solving the equations of motion 

Let us now try to solve the equations (|6.1 |) - (|6.5|) in Fourier space. We propose 



A\x) 
N\x) 
L\x) 



fAkT. {AU'l^t)e(-^'i<r)e"^') +h.c., 

^ ' m=l 



■J \ J m=l 



b{x) = I ^[-doMQ,t)-i\q\Mq,t))e^'^^ + h.c.. (6.18) 

The next step is to derive the equations for the Fourier coefficients. Plugging (|6.18|l into (|6.2|l , 
we derive 

{di+f)Ao = mi\q]N3, (6.19) 

while ( |6.3I > yields 

{dl + f + m^)Ns = m (doAs + i\q\Ao^ , (6.20) 
(di + f + m^)Ni,2 = m5oIi,2. (6.21) 

(|6.4|l results in 

3 3 

-{di + f + m'^)J2^rn£^"'^'' = me^^' V' ^ e^"^)'*^!™ . (6.22) 

m=l 771=1 

Multiplying with e^")-* , summing over i and using the orthonormality, we find 

3 

- {d^ + <f + m')Ln = -im Amq' (^^"^ x e^^^y . (6.23) 

m=l 

Taking n = 2>, we know that will only survive for m = 1, 2, however then the cross 

product is oc e(2)'(^) and thus orthogonal to q, meaning that 

{dl + q^ + m'^)L^ = 0. (6.24) 

With n = 1, 2, dESll) is valid if 

- {dl + q^ + m^)Li = -im\^A2, (6.25) 
-{dl + q^ + m'^)L2 = im\q\M . (6.26) 
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Analogous manipulations on (|6.5|) give 

{dl+q^)A^ + mdoN^ = 0, (6.27) 

{dl + q^)Ai + mdoNi = -im\q\L2, (6.28) 

{di+q'^)A2+mdoN2 = im\^Li . (6.29) 

Let us now try to find a sensible solution to the previous differential equations in t. 

We first concentrate on the independent subset of equations (|6.19|l , (|6.20|) and (|6.27|l . We 
can decouple these differential equations by passing to higher order differential equations. 
These induced equations read 

{di + ff{di + f + M^)Ao,3 = 0, 

{di + (f){di + f + M^)Ns = 0, (6.30) 

where we defined 

= + . (6.31) 

We solve these equations bjS 



mq,t) = n3(g)e-*l'1*+n3(g1e-V^^*. (6.32) 

For consistency, we must impose the original set ( I6.19D , (I6.20D and ( I6.27D again, and identify 

the coefficients of the various t-dependent functions (e~*l^*, te"*''^* and e^*V^^^+^*). The 
following equations come out as independent ones 

— M^ao = mi\q\ns , 

—2aQ = mns , 

as = oo , 

—mi\q]as + mas + milqlat) = m^ns , 



-M^ag = miy/M^ + q^n^ , (6.33) 

implying that there are only 3 independent coefficients left of the original 8 in (|6.32|l . 
The next step is to analyze the remaining equations. (|6.24|) immediately gives 

Mq, t) = i3{q)e-'^/^+^' . (6.34) 

The last independent set consists of (|6.2H) , (|6.25|) , (|6.26|) , (|6.28|) & (|6.29|) . Applying the same 
trick as before, we deduce 

{di + f){di + f + M^)A,^2 = 0, 
{di + (f){di+f + m^){di+f + M^)N,,2 = 0, 
{di + f){di + f + m^){di + f + M^)L,,2 = 0, (6.35) 



We did not write explicitly the h.c. part. 
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with corresponding solution^ 



L2(g, t) = h{q)e-'\'^' + ^{q)^-'^^^'' + h{q)e-'^^^' , 

N2{q.t) = n2{q)e-'\'^' + n2{q)e-'^^^' + n2{q)e-'^^^' . (6.36) 

Consistency requires 

fh?ni = —im\<f\ai , fh?n2 = —im\<f\a2 , 

mni = i\J (p' + M^ai , mn2 = i\J (p + M'^a2 , 

mil = —'i'\q\a2 , rn(^2 = ^l^loi , (6.37) 

ni = i2 , _ n2 = -h , 

This reduces the number of independent Fourier coefficients in (|6.36|) from 16 to 6. 

Summarizing, we have 3+1+6=10 independent Fourier coefficients left (plus their hermitian 
conjugates). Without loss of generality, we choose to work with n3{q), ai{q), a2{q), ai{q), 
«2(g), n3{q), £3(0), ii{q), £2(0) and ao{q). 



6.2 The conjugate momenta 

In order to quantize the theory, we need the conjugate momenta. We shall only be concerned 
about the degrees of freedom hidden in the bosonic fields. We do not care about the ghosts 
for the moment. 

Making use of the action (15.121 ), we derive the desired conjugate momenta. 



6C 



< = T77n^ = -K-mNt, 



7r°'" = \f- = -d^Al. (6.38) 
For later use, let us compute the the multiplier h of (|6.18|) , which is actually given by 

m2 r 



b{x) 

ni 



I (^"3(g1e-^l«1*e^'?x+h.c.. (6.39) 



^We assume that / 0, otherwise we should adapt the analysis. 
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Let us now quantize the model. Naively, the brackets we would impose arj^ 

[Ai{x,t),Trj{y,t)] = -i6ij6^^\x - y) , 
[Ao{x,t),7T'>{y,t)] = i6('Hx-y), 

|c(x,t),7rc(jr,i)} = i6^^\x-y), 

|c(f,t),7rc(y,t)} = i6^^\x-y), 
[N,{x,t),7rf{y,t)] = -idij6^^^x - y) , 

[Li{x,t),7rf{y,t)] = i6,,6'^^\x - y) . (6.40) 

However, these brackets will be not necessarily correct, due to the fact that multiple relations 
exist between the different field and conjugate momenta configurations, something which 
is clearly visible from the relations between the Fourier coefficients (|6.33l > and (16.371 1. We 
postpone the actual discussion of the quantization to section 7, where we shall explain an 
alternative way to fix the commutation relations in an appropriate fashion. 

An additional complication that we should have a look at is the precise form of the free 
Hamiltonian Ho- In order to have a well defined Fock state space describing physical parti- 
cles, the states we are considering, which are defined by acting with the creation operators 
on the vacuum, should be energy eigenstates of the Hamiltonian, which should therefore be 
diagonalized in terms of certain creation and annihilation operators. 



6.3 Calculation of the free Hamiltonian 

We shall neglect that part of the Hamiltonian depending on the Faddeev-Popov ghost fields 
7f and c", as well as depending on the fields a*'°, a*'"^, x'^''^'^ and these are 

of no relevance for the present discussion, and we call this "reduced" Hamiltonian Hq. Since 

K = J d"^^ [nidoA' + 690^0 + vr^,i5oiV' + 7rL,idoL' - C] , (6.41) 



we find 



n'o = j d^x 



-bdf.A^' - + + mN')doA' + bdoAo + (^oA^* - mA')doN' 



- doL'doV - ^F'^^F^' + ^F'^F'^ - mN\d^A' - d'A^) + me'^^L'^d'A^ 

- -doN'doN' + -N'djd^N' + -doVdoV - -Vdjd^V - -m^VV + -m'^N'N' 
2 2 2 2 2 2 

(6.42) 

Since we must plug in the solutions (|6.18|) into the field expression of the Hamiltonian Hq to 
find its operator valued expression, we can already use the equations of motions ( I6.1| )-( I6.5D 



^^See also the Appendix. 
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to simplify a bit H'q. Since 



F 



Oi 



d^A' - d'A 



and 



ppg _ Qp^g _ gg^p 
=^ {dPAi - d^APf 



m 
1 

m 
1 



+ m^) N' , 
+ m2) Lk 



m 



(6.43) 



we simplify H'q to 



2 m 



- doVdoL 

- -doN'doN' + -N'djd^N' + -doVdoL' - -Vdjd^V - -m'^VV + -m'^N'N' 
2 2 2 2 2 2 

(6.44) 

As a first exercise, let us have a look at the terms that give rise to the ^s-oscillator. Due to 
(|6.24|) , it is easy to see that the only terms relevant for this are 



d' 



9 



{(f + "1.^)^3 (5)^3 (g) after normal ordering . 



(2vr)^ 



(6.45) 



We immediately notice that we have a negative sign in the Hamiltonian as far as the £3- 
modes are concerned, which is inconsistent with a positive commutator for the reason of 
positivity (see also section 7). Invoking a negative commutator then necessarily leads to 
negative norm states. It might be important to notice here that the BRST charge (|4.5|l is of 
no help whatsoever to eliminate the negative norm states created by the oscillator ^. The 
charge Qq does only act on massless states, as implied by (|4.5|l and l|6.39|) . 

One might wonder if there might exist a way out, in order to find a positive Hamiltonian, 
without the need for negative norm states. If we could change the sign of the "new" terms 
in the gauge model, we should at least be able to avoid the problem in the £s sector. 
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We would thus like to use the following action 

Sd + Sgf , 



^phys 



(B - B) F , ^ ^ ^ ^ 



(6.46) 

ji'^ T^ab T^bc Tjc j^ab j^bc^c \ 



' 111/-' 



32 



+ 



\abcd 



16 



B^Bt, - GLct,] ( Bint - Gl^d 



fiv flU ^ fii/^ fiv J \ pa pa ^ pa^ pa 



(6.47) 



which is completely similar to the action up to a few signs. 

The steps in the algebraic renormalizability analysis of 1TH2| are not essentially affected by 
these sign changes. More precisely, it is still possible to show the renormalizability to all or- 
ders. It is clear that the supersymmetry ( |5.6I > is still present for m = 0, so that the equivalence 
with ordinary Yang-Mills theories is maintained for m = 0. 

In the decomposed form, we have the following (quadratic) action in Minkowski space time 



So 



= J d'^x {df,Al - d^Alf - mN'''' {d^A'" - d'A^") + me'^^L^^" {d'A^") 
2 2 2 2 2 

- rr?^ (Ar*'«Ari.« _ Li,a^i,a-^ ^ b'^d^A'"' + ^d'^c" + ^b^b" . (6.48) 



Evidently, the cohomology analysis of section 4 can be immediately translated into the new 
language. 

The classical equations of motion read 



b + df'A^ = 0, 
d'^Ao + mdiN' = 0, 



0, 



-{d^ + m'^)L' -me'^'^djA'' = 0, 
d'^A^ -mdoN^ + me'^^djL^ = 0. 



(6.49) 
(6.50) 
(6.51) 
(6.52) 
(6.53) 



We solve them by using once more the Fourier decomposition (|6.18|) , yielding the following 
differential equations in time t 



{dl + t + rr?)N-i = m[dQA^ + i\q\AQ 
{dl + q^)A^ - mdoN^ = 0, 



(6.54) 
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and 



{dl + f + m'')Ni^2 = mdoAi^2, 

{dl+q^+m^)L^ = 0, 

(^0 + m^)Li = im\q\A2 , 

{dl + + ?iT?)L2 = -im\q\Ai, 

{dl + q^)Ai - mdoNi = im\^L2, 

{dl + q^)A2 - mdoN2 = -im\q\Li . (6.55) 



Decoupling the first set of equations (|6.54|) leads to 



{dl + ff{dl + t + M^)A^,^ = 0, 

{dl + if){dl + f + M^)N^ = 0, (6.56) 

where the mass is now defined by 

The presence of the mass scale m? 7^ is essential here to avoid the appearance of a tachyonic 
mass —m?, which is due to the change of sign in the action (|6.46|) . The mass scales oc Xim? 
leading to the mass m? were originally introduced in [IJ to ensure that the final action is 
renormalizable from the algebraic point of view. In the present context, these also play a 
physical role in order to avoid tachyons, if we assume that m? > m?. 

The solution to (|6.56|) is still given by (|6.32|) , however the relations between the coefficients 
are modified into 

M^oo = mi\q\n^ , 
2ao = mn3 , 

03 = So , 
—mi\q\a^ + ma-^ + 'mi\q\aQ = in^n^ , 



M^as = miy^M^ + g^ns . (6.58) 



It is straightforward to verify that 



bix) = /(^(^)"3(g>-^I^^V-+h.c.. (6.59) 
The conjugate momenta are given by 
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The equation of motion for L3 is immediately solved by (|6.34|) . 

Let us now turn to the (1, 2) -sector. It is not difficult to check that (|6.35|) and (|6.36|) still holds. 



but now with 



m?n2 = — im|g|a2 



m?ni = —im\q\ai , 

mni = + M^ai , mn2 = —i\/<p + M'^a2 , 

mil = ^1^02 , m£2 = -i\q\ai , 

ni = £2, n2 = -h , 



(6.61) 



The (reduced) free Hamiltonian H'q becomes 
Ho 



2 4 ^ ^ 

2 2 2 2 2 2 



1 



(6.62) 



Since the relations (|6.43|) are still valid, we find for the first part 



I on-shell 



-b^ + — (9^ + M'^)N^dQA^ + 65o^o + mA'd^N^ 
2 m 



2m^ 2m^ 



1 



N'did^N' 



(6.63) 



',(1) 



After some tedious algebra, we can express the part of the Hamiltonian corresponding to the 
(massive) '^-modes as 

d^q 
(2^ 

- 2{f+ m^) n{{q)ni{q) - 2 {f + m^) n\{q)n2{q) 
Using the relations {6.5%} and (|6.61|) , we rewrite (|6.64|) as 



(6.64) 



',(1) 



d\ r, 

(27r)= 



2 (g2 + ^ (^^-^ ^ 2m2^T (g)4 (g) + 2?f?Il {q)£2 {q) ■ (6.65) 



2^/ 



The second part, corresponding the (massive) A-modes is given by 



',(2) 



d^q 



1 \ 



2(g^ + M^) ll + j-^ \ nl{q)n^{q) - 2m^n\{q)ni{q) - 2vfi'n\{q)n2{q) 



;2c.t/ 



) \ 



(6.66) 



Finally, we also have a part 'H^^^ corresponding to the massless modes, which we do not 
write in more detail. 
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7 Fixing the commutation relations between creation and annihi- 
lation operators and the non-unitarity of the model 



7.1 Outline of the idea 

We recall here that the equations of motion (|6.49|) - (|6.53|) constitute several (nontrivial) rela- 
tions between the field components. Unfortunately, these relations are quite complicated 
to solve explicitly. Closely related to this is the fact that it is not immediately clear how to 
diagonalize the quadratic form appearing in the action (|6.48|) . Moreover, this would-be di- 
agonalized form should also be of maximum second order in the derivatives to allow for a 
consistent quantization. 

Since the fields are related to each other, we can expect the same for the conjugate momenta. 
E.g., it is easily checked that vro(x, t) = —7r^{x, t). 

It would appear that quantization, starting from configuration space, is a highly nontrivial 
task due to the existing relationships between the fields and momenta. It would be beneficial 
to quantize the theory directly in Fourier (momentum) space, in which case the independent 
degrees of freedom are clearly identifiable. 

There does exist a simple way to derive the appropriate (anti-)commutation relations be- 
tween the creation and annihilation operators, when we rely on the Heisenberg picture. We 
can always calculate the free Hamiltonian operator Ho without exact knowledge of the com- 
mutation relations. The precise value of the commutators only influences the constant part 
of Ho, which we can drop after normal ordering. On one hand, we already know the time 
evolution of all fields as written down in (|6.18|) by using the Fourier decomposition and 
the solutions (|6.32|) and (|6.36|) . On the other hand, these fields are treated as operators in 
the Heisenberg picture, and consequently their time evolution is dictated by the Heisenberg 
equation 



Demanding consistency allows to fix the commutators between the creation and armihilation 
operators, which are supposed to be number valued. As an illustration, consider the Klein- 
Gordon action 



[Ho,^] = -i^^^. 



(7.1) 




(7.2) 



The solution to the equation of motion is expressed by 




(J, 



(7.3) 



As it is well known, the normal ordered Hamiltonian is given 



b30 




(7.4) 



14- 



'We disregard the conventional : Hkg '■ notation and continue to use Hkg- 
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so that requiring 



d 

[Hkg, ^{x, t)] = -i—ip(x, t) , 



leads to the correct commutation relation 



(7.5) 



(7.6) 



7.2 Quantization of the massive gauge model using the Hamiltonian TYq 

In order to find the commutators for the massive '^-modes, we only need to consider T-Iq^^\ 
When we impose 



we deduce that the following commutation relations arise 

X,2(g),?i,2(^) 



{2nf 
(2vr)3 



2m? 
1 



2^/W+ 



others trivial 
We can define novel operators as follows 



2\J (p + rri 



to arrive at a free Hamiltonian in a more standard form 



',(1) 



with commutation relations 



d?q 1 
(27r)3 25^ 



2^,Y.^'l{q)m 



i=l 



others trivial . 
The energy was defined as 



Zjq{2T:f5ij5^^\k 



(7.7) 



(7.8) 



(7.9) 



(7.10) 



(7.11) 



(7.12) 



We conclude that the modes £'^{q) correspond to the 3 polarizations of a vector particle with 
mass m. As it is clear from (|7.1H) , they have a positive norm. 

If we would have used the original Hamiltonian (|6.44|) , we would have found a negative 
sign commutator, as already advocated below (|6.63|) . 
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In a completely similar fashion, we can quantize the A-modes. We rewrite (|6.66|) into the 
form 



n 



',(2) 



d^q 1 
(27r)3 22q 



i=l 



(7.13) 



by a suitable rescaling of the fields, in which case the corresponding commutation relations 
become 



others trivial . 



with 



(7.14) 



(7.15) 



We conclude that we have been able to change the negative norms corresponding to the 
£i-modes in our first model (|1.1|| by passing to (|6.47|) , but unfortunately the negative norm 
problem has been "shifted" to the rij-modes, as it is clearly visible from (|7.14|l . We are thus 
forced to conclude that our theory is not unitary, as we do not seem to be able to remove these 
negative norm states, corresponding to the polarizations of a massive vector particle with 
mass M. Just as in the case of the ^j-modes, we reckon once more that the BRST symmetry 
has nothing to say about the fij-modes. 

To completely finish the analysis, we should in principle also determine the free Hamilto- 
nian corresponding to the massless modes. However, this has become a bit redundant job by 
now since we already have established that the theory cannot be unitary. Quite obviously, 
some restriction on the allowed (massless) states and Faddeev-Popov ghost states will arise 
from the BRST charge (|4.5|) . However, this might also turn out to be a not so trivial task after 
all due to the multipole character, encoded in (|6.32|) in the terms linear in time t. As it was 
found in the case of the toy model of L46J, the Hamiltonian is not even necessarily diagonal, 
obscuring the particle interpretation. It is even noticed that the Hamiltonian cannot be diag- 
onalized with linear transformations. As already said, we shall however not dwell on that 
point here. 



8 Conclusion 

We continued our investigation of a recently proposed local non-Abelian gauge invariant 
action containing mass terms il.l} , which can be renormalized to all orders of perturbation 
once a gauge is chosen. Despite the fact that the model enjoys a BRST symmetry with nilpo- 
tent generator given in (|1.14|) , it turned out to be impossible to remove all negative norm 
states from the physical state space. As a consequence, our model is not unitary and thus 
not useful as a physical theory containing massive particles as asymtotic observables. 

As already indicated in the introduction, this does not mean that our action is now useless. 
We can still wonder whether it could be generated dynamically in a nonperturbative fashion 
when we would start with the initially massless version of our gauge model. As a matter 
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of fact, we can propose our model as an alternative to ordinary Yang-Mills theory in the 
perturbative (massless) region, since both are equivalent [2J. It might however be more 
convenient to search for a dynamical mass m in our model since the mass term can be at least 
coupled to the action without spoiling the renormalizability or gauge invariance. Of course, 
it remains to be investigated if a sensible gap equation could be established. If this would 
work out, we might have a gauge invariant mechanism behind the dynamical generation of 
a massive parameter into a gauge theory. 

We repeat that the generation of a mass would not be necessarily in conflict with the nonuni- 
tarity of the massive gauge model, since we start from the (unitary) massless theory and we 
do no longer want to describe the asymptotic high energy behaviour wherefore perturbation 
theory applies perfectly, but rather we are entering a phenomenologically interesting region 
where e.g. the gluons already loose their physical meaning as an observable. 

Since our generated mass would be gauge invariant, it could enter physical correlation func- 
tions. As such, it could be investigated if it could serve as a possible alternative to the {A'^)mm 
condensate used to explain the ^ power corrections IBIZl- 

We conclude by noticing that we could repeat our analysis in the Abelian case. Then it can 
be shown that the action il.l^ is stable without the need for Ai and A3 UJ. It can also be 
shown that the model is equivalent with the Abelian Stueckelberg model ||47|. described by 



When the auxiliary fields are integrated out in both cases, some manipulation leads to the 
same (nonlocal) action 



It is known that the Abelian Stueckelberg model is renormalizable and unitary, see e.g. l47| 
for a review. If we would analyze the unitarity of the Abelian version of (|1.2|), we would run 
into exactly the same problem as in the non- Abelian case, i.e. the presence of negative norm 
states in the physical subspace. 

The lesson to be learnt is the following. We depart from the same nonlocal action, but in 
order to give a consistent quantization of the theory, including an analysis of the unitarity 
and renormalizability, we are forced to bring the action in some localized polynomial form. 
Apparently, the exact procedure of localization affects these results. For example, in the 
Abelian case, there are 2 quite distinct approaches to bring the action (|8.2|) in a local and 
in addition renormalizable form: the Stueckelberg approach or the pathway we followed. 
However, only the Stueckelberg way gives a unitary model. Apparently, the precise role of 
the additional fields that are introduced cannot be underestimated in the discussion of the 
unitarity and/ or renormalizability. 




(8.1) 



'nonlocal 
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Appendix 

Since we have been working with the decomposed fields, the covariance of the (naive) com- 
mutation relations (|6.40|) is obscured. It is worth having a look at this. We decompose 
and its complex conjugate B^^ into their real and imaginary part. 



Xf_iu = -^{B^u + B^u) Y^u = '^{Bfiu — B^i,) , (8.3) 



then 



Mi = Xoi , Ni = Yoi , 

= ^e'^'^Xjk , = \e'^%k ■ 

Covariance and the antisymmetry would require commutation relations like 

[X^i.(x,t),7r^^(y,t)] = i{giicc9v(i- gii(igva)f>^^\x -y) , 



Specifically 



(8.4) 



[y^i.(x, t), 7r^^(y, t)\ = i {g^,c9uf3 - gfifsgua) S^^Hx - y) . (8.5) 



[Xoi(x,t),7r^(y,t)] = igijd'^^^x - y) = -i6ij5^^\x - y) , 
[Xij(f,t),7r^(y,t)] = i{5^k5jl-5il5Jk)5^'^Hx-y), (8.6) 

and the same for Y . Consequently 

[M,(f,t),7rf (y,t)] = -i6ij5^^\x - y) , 
[Kr{x,t),TTf{y,t)\ = i£"%^'='[X,,(x,t),7r^(y,t)] 

= i^e"%^'=' {5ik5,i - 5u5,k) S^^^ {x-y) = i6rsd^^^ {x - y) , (8.7) 

and the same for N and L. 
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